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The thermodynamics of the SU(3) gauge theory has been analyzed with tree level and tadpole improved 
Symanzik actions. A comparison with the continuum extrapolated results for the standard Wilson action shows 
that improved actions lead to a drastic reduction of finite cut-off effects already on lattices with temporal extent 
N T = 4. Results for the pressure, the critical temperature, surface tension and latent heat are presented. First 
results for the thermodynamics of four-flavour QCD with an improved staggered action are also presented. They 
indicate similarly large improvement factors for bulk thermodynamics. 



1. Why improved actions are expected to 
improve the QCD thermodynamics 

Although the plasma phase of QCD does 
in many respects show strong non-perturbative 
properties (screening masses, poor convergence 
of perturbation theory,...), bulk thermodynamic 
observables like energy density and pressure do 
approach the non-interacting ideal gas limit at 
high temperatures (T) where they are expected 
to receive large contributions from high momen- 
tum modes, (avg. momentum) ~ T. In lattice cal- 
culations T is determined through the temporal 
extent N T of the lattice and the cut-off a -1 , i.e. 
T = l/N T a. In the ideal gas limit the relevant 
momenta are therefore of the order of the cut-off, 
where lattice and continuum dispersion relations 
differ strongly from each other. Indeed this is the 
origin of the well-known discrepancy between the 
energy density of an ideal gas calculated on a fi- 
nite lattice (e(N T )) and in the continuum (csb), 
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The standard Wilson action leads to nearly 50% 
corrections on a lattice with temporal extent 
N T = 4 while the cut-off dependence is drastically 
reduced in the case of Symanzik improved actions 
which eliminate the leading 0(N~ 2 ) = 0{(aT) 2 ) 
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corrections. In the case of the (1 ,2)-Symanzik 
action this is less than 2% for N T = 4, i.e. 
ci = 0.044. This clearly demonstrates the signif- 
icance of improved actions in the ideal gas limit 
(T-oo). 

In a first exploratory analysis we have shown 
that the improvement found analytically in the 
ideal gas limit also persists at finite temperature 
and does seem to be important even close to T c 
0. Here we will present further evidence that 
tree level and tadpole improved actions do lead 
to a significant reduction of finite cut-off effects 
even at T c §-g. 

2. SU(3) Thermodynamics 

We have analyzed the thermodynamics of the 
SU(3) gauge theory using 0(a 2 ) improved tree 
level and tadpole improved actions, which have 
been constructed by adding to the Wilson pla- 
quette action, S^ 1 ' 1 ), an appropriately weighted 
contribution from planar (1,2) or (2,2) loops, 
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Table 1 

Critical temperature in units of y/a on lattices with temporal extent N T — 4. Infinite volume extrapola- 
tions of P c for the Symanzik actions are based on lattices with N a — 16, 24 and 32. 



action 


N a 


0c 




standard Wilson 


oo 


5.69254 (24) 


0.5983 (30) 


(2,2) Symanzik (tree level) 


24 


4.3995 (2) 


0.624 (4) 


(1,2) Symanzik (tree level) 


oo 


4.07297 (28) 


0.631 (3) 


(1,2) Symanzik (tadpole) 


oo 


4.35228 (39) 


0.635 (3) 



where the factor uq appearing in the definition 
of the (1,2)-Symanzik action denotes the tad- 
pole improvement factor defined in terms of the 
self-consistently determined plaquette expecta- 
tion value @] . 

2.1. Bulk thermodynamics 

Bulk thermodynamic quantities like the energy 
density (e) or pressure (p) calculated with the 
Wilson action on lattices of size TV 3 x N T with 
N T = 4, 6 and 8 || show a strong cut-off de- 
pendence in the plasma phase. A first analy- 
sis with the (2,2)-Symanzik action jl| has shown 
that this cut-off dependence gets drastically re- 
duced already on a N T = 4 lattice. The work 
presented in Ref. |jj has been extended in several 
ways. On lattices with temporal extent iV r = 4 
we have analyzed now also the tree-level im- 
proved (1,2)-Symanzik action, which shows even 
smaller cut-off distortions in the high tempera- 
ture limit than the (2,2) action. In addition we 
have studied the influence of tadpole improve- 
ment in the case of the (1,2)-Symanzik action. 
In all cases we have determined the temperature 
scale non-perturbatively through the calculation 
of the string tension on symmetric (N*) lattices, 
T/T c = y/aa(f3 c )/ \faa(0). From this calculation 
we also obtain T c in units of y/a for the different 
actions. These are given in Table [j] together with 
results for the Wilson action. In that case T c /y/a 
has also been calculated on lattices with N T vary- 
ing between 4 and 12 and extrapolated to the con- 
tinuum limit. This gave T c j 'y/a = 0.629 ± 0.003 
|| . On iV r = 4 lattices the results obtained with 
improved actions are thus much closer to this ex- 
trapolated value than those obtained with the 
Wilson action. 



The pressure has been calculated by integrating 
differences of the action densities obtained on 16 4 
and 16 3 x 4 lattices (in the case of the (2,2) action 
we used 24 4 and 24 3 x 4 lattices) following the 
standard procedures Jl]||. Results are shown in 
Figure [j]. 
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Figure 1. Pressure of the SU(3) gauge theory on 
lattices with temporal extent N T — 4. The solid 
line shows the continuum extrapolation obtained 
from the standard Wilson action. The dots give 
results from a calculation with a perfect action 
on a 12 3 x 3 lattice [6]. The arrow indicates the 
continuum ideal gas value. 

2.2. Surface tension and latent heat 

The success of improved actions for the calcu- 
lation of bulk thermodynamics even at tempera- 
tures close to T c naturally leads to the question 
whether these actions also do lead to an improve- 
ment at T c . This is, of course, a highly non- 
perturbative regime. However, observables like 
the latent heat (Ae) and the surface tension (a/), 
which characterize the discontinuities at the first 
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order deconfinement phase transition in a SU(3) 
gauge theory, do depend on properties of the low 
as well as the high temperature phase. As the lat- 
ter is largely controlled by high momentum modes 
it may be expected that some improvement does 
result even from tree level improved actions. 

At the SU(3) deconfinement transition Ae and 
07 have been studied on lattices up to temporal 
extent N T = 6 |||| . A strong cut-off dependence 
has been found when comparing calculations for 
N T = 4 and 6. For this reason an extrapolation to 
the continuum limit has so far not been possible 
for these observables. 

We have extracted erj || from the probability 
distribution of the absolute value of the Polyakov 
loop following the analysis presented in Ref. ||. 
The probability distribution at the minimum is 
proportional to 

P(\L\) ~ exp(- + f 2 V 2 + la 1 A] /T) (3) 

where /j denotes the free energy in the phase i, 
and Vi is the volume occupied by that phase and 
A denotes the interface area of the finite system. 
From the depth of the minimum one thus can ex- 
tract the surface tension. The distribution func- 
tions for the tadpole improved actions for three 
different lattice sizes are shown in Figure ^. 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 

Figure 2. Polyakov loop distributions for the tad- 
pole improved (1,2) Symanzik action. 

In Table || we give results for aj on the largest 
lattices considered. Clearly the surface tension 



extracted from simulations with improved actions 
on lattices with temporal extent N T = 4 are sub- 
stantially smaller than corresponding results for 
the Wilson action. In fact, they are compatible 
with the N T = 6 results for the Wilson action. 

The latent heat is calculated from the discon- 
tinuity in (e — 3P). This in turn is obtained from 
the discontinuity in the various Wilson loops en- 
tering the definition of the improved actions, 

I - ;(&)'(•£)(<*-<*-)•<" 

with S = S- &S/d(3. 

The difference of action expectation values at 
f3 c is obtained by calculating these in the two co- 
existing phases at (3 C . In the time histories of the 
Polyakov loop values we have cuted out the tran- 
sition regions in order to classify configurations 
belonging to either of the two phases || . 

In order to extract the latent heat one does 
still need the /3-function entering the definition 
of Ae/T^ in Eq. || The necessary relation a(/3) 
has been obtained from a calculation of \faa (im- 
proved actions) or a determination of T c a (Wilson 
action). Results for Ae/T* are summarized in Ta- 
ble H Here we also give the result obtained with 
the 1-loop /3-function. 

3. Four-flavour QCD with an improved 
staggered action 

In the fermionic sector of QCD the influence 
of a finite cut-off on bulk thermodynamic observ- 
ables is known to be even larger than in the pure 
gauge sector. For instance, in the staggered for- 
mulation the energy density of an ideal fcrmi gas 
differs by more than 70% from the continuum 
value on a lattice with temporal extent N T = 4 
and approaches the continuum value only very 
slowly with increasing N T . This cut-off depen- 
dence can drastically be reduced with an 0(a 2 ) 
improved staggered action. We use an improved 
action, S'lU] = S^ 1 ' 2 ) + 'tjjMijj, where a higher 
order difference scheme (one-link and three-link 
terms) , is used to improve the fermionic part Jl0[ . 
The improved fermion matrix reads 

Mp^ = mS id + m ( ^-A [U]. j - ±-B [C7]. .) 
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Table 2 

Surface tension and latent heat for the three improved actions and the Wilson action. Results for the 
Wilson action are based on data from [9] using the non-perturbative (3- function calculated in [5]. 



action 


v a 


N T 




Ae/T c 4 


(a £ /t c 4 ) 

/ pert 


standard Wilson 


24 2 x 36 


4 


0.0300 (16) 


2.27 (5) 


4.06 (8) 




36 2 x 48 


6 


0.0164 (26) 


1.53 (4) 


2.39 (6) 


(1,2) Symanzik (tree level) 


32 M 


4 


0.0116 (23) 


1.57 (12) 


2.28 (8) 


(1,2) Symanzik (tadpole) 


32 3 


4 


0.0125 (17) 


1.40 (9) 


1.94 (8) 



with 

J" 

B[U\ij = ^^(CA iAt {7j_|_p, )At J7 ! i+2A,Ai^ij-3A — 

~^i-(i,pPi-2jl,pPi-?,[L,^ i ,i + 'iiJ-) 

With this action the overall cut-off distortion of 
the ideal gas limit on a 16' ! x 4 lattices reduces to 
about 20%. We have performed simulations for 
two quark masses, ma = 0.05 and 0.1 Like in 
the pure gauge sector the improvement is visible 
already close to T c . For instance, the energy den- 
sity stays close to the ideal gas limit. Although 
we observe an overshooting of the ideal gas limit 
close to T c for the non-zero quark masses con- 
sidered by us this feature does seem to disappear 
in the "chiral limit" {§. 
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Figure 3. Energy density of four-flavour QCD on 
a 16 3 x 4 lattice. The lower set of curves shows 
an extrapolation to the "chiral limit". 



4. Conclusions 

Thermodynamic observables of SU(3) gauge 
theory and QCD studied with improved gauge 
and fermion actions show a drastic reduction of 
the cut-off dependence in the high temperature 
limit as well as at T c . The major improvement 
effect is already obtained with tree level improved 
actions. Tadpole improvement plays a minor role 
even close to T c . 
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